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Particle Loss in Sampling Lines 

Thomas Wey and Nan-Suey Liu 
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Glenn Research Center 
Cleveland, Ohio 44135 


Summary 

This paper at first describes the fluid network approach recently implemented into the National 
Combustion Code (NCC) for the simulation of transport of aerosols (volatile particles and soot) in the 
particulate sampling systems. This network-based approach complements the other two approaches 
already in the NCC, namely, the lower-order temporal approach and the CFD-based approach. The 
accuracy and the computational costs of these three approaches are then investigated in terms of their 
application to the prediction of particle losses through sample transmission and distribution lines. Their 
predictive capabilities are assessed by comparing the computed results with the experimental data. The 
present work will help establish standard methodologies for measuring the size and concentration of 
particles in high-temperature, high-velocity jet engine exhaust. Furthermore, the present work also 
represents the first step of a long term effort of validating physics-based tools for the prediction of aircraft 
particulate emissions. 


Introduction 

Aircraft emissions contribute to climate change in the atmosphere and air pollution in the 
neighborhood of airports. Assessment of the magnitude of their environmental impacts calls for a better 
understanding of the formation and the subsequent development of gaseous pollutants, aerosols (volatile 
particles and soot), and their precursors in the internal and in the plume of the jet engines operating over 
the full range from ground to flight altitude. The volatile nature of aerosols introduces difficulties in their 
repeatable and accurate measurement because their concentration and size can be quite sensitive to the 
sampling conditions and procedures. These difficulties are sometimes further compounded by the 
necessity of in-situ measurements, i.e., sampling from aircraft in flight. Clearly, the ability to predict the 
effects of sampling conditions and procedures on the aerosol behavior will greatly assist the development 
of sampling techniques and the optimization of sampling procedures for making quantitatively accurate 
measurements of particulate emissions. 

A multi-dimensional CFD-based capability for modeling and simulation of the aviation-sourced 
particle evolution in flow paths is now available in the National Combustion Code (NCC). A detailed 
description of this multi-dimensional approach can be found in references 1 and 2. The CFD-based 
simulation is very computing intensive, and, quite often, a lower-order temporal approach is used to 
capture the overall evolution, a detailed description of this approach is given in reference 3. This lower- 
order temporal approach greatly alleviates the computing resource requirement, but it may not be 
appropriate for certain flow configurations and conditions. To bridge the gap between the CFD-based 
approach and the lower-order temporal approach, we have implemented the fluid network approach into 
the NCC, the details of the adopted fluid network approach will be presented in the next section. 

The NASA-sponsored Aircraft Particle Emission Experiment (APEX) field test series to measure the 
size and concentration of particles in the jet engine exhausts have produced highly variable results (ref. 4). 
To better understand the sources of this variation, NASA recently sponsored additional field experiments 
and laboratory experiments to evaluate jet engine exhaust sampling and measurement systems. One of the 
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focused area is the assessment of size and mass-dependent particle loss through the sampling lines. These 
tests produced several sets of particle line loss data, which are used in the present effort to validate the 
predicted results obtained from the three different modeling and simulation approaches mentioned above. 
Comparison between the measured data and computed results are discussed in a later section. 


Modeling and Simulation 

Same Chemical and microphysical models for aircraft particulates are used in the three different 
levels of approach. The finite rate kinetics models used are those for conditions intermediate between 
combustion and atmospheric chemistry. The particle transport equations are based on the multi-bin 
Eulerian framework. The microphysical models include sulfate and nitrate aerosol precursors, 
thermophoresis, coagulation, nucleation, soot activation, condensation/evaporation, and interaction 
between the gaseous species and the particles. Generally speaking, the CFD-based simulation tightly 
integrates fluid dynamics, chemistry, and particle microphysics. The lower-order temporal simulation is 
one in which chemistry and particle microphysics are driven by averaged flow parameters specified as a 
function of time. Details of these two approaches have been reported in references 1, 2, and 3. 

Motivated by the practical need to simulate the often complicated sampling line system in an efficient 
and timely manner, we added an additional option based on the fluid network approach, which takes into 
account the complex internal flow and thermal effects within the fluid systems. Although here the 
intended application is for the particle sampling line system, the fluid network approach is, in general, a 
system level method and is not restricted to any specific geometries or configurations. SINDA/FLUINT 
(ref. 5) and GFSSP (ref. 6) are two of many system simulation programs widely used by NASA for 
internal flows. This work adopts similar formulation of the momentum equation as that in GFSSP. A 
relatively detailed description of the fluid network approach used in this work follows. 

The conservation equations for the gaseous flow and the transport equations for the aerosols are 


Continuity equation 

iSISr pdV+ S 1^-^=° 


Species transport equation 

dV + \\ A PmU ' dA = \\ A pDe ff V(R p ' dA + JJJ, vv ">» dV 


( 1 ) 


(2) 


Momentum equation 


\\1 PU dd + \\ 4 Pdd ' ~ ~ || ( pdA + 1| ( r ■ dA 


( 3 ) 


Energy equation 

\\\v PE dE + \\ 4 PEd ' dA ~ ~ \\ A pd dA ~ || | ^ ^ + \\ A U ' T 'dA (4) 

Aerosol general dynamic equation 

flfl> dV + \\ A N k" ■ dA = \\ A iD M + -jh=)VN k ■ dA 
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where 



( 5 ) 


V is an arbitrary control volume with control surface A, p is the fluid density, p„, is the partial density of 
species m, \ Vm is the production rate of species m per unit volume, u is the flow velocity in a stationary 
Cartesian coordinate system, E is the total energy density, p is the static pressure, t is the viscous stress 
tensor, q is the heat flux vector, D e g\s the effective mass diffusivity coefficient which is equal to ' U ^‘“ , 

k, p, and p, are the thermal conductivity, laminar viscosity and turbulent viscosity, respectively. For gas, 
the Schmidt number (Sc) is about 1. Nk is the number density of particles in size bin k, M B is the total 
number of the size bins, v, is turbulent kinematic viscosity, C u is 0.09,/, is the damping function of the 
turbulence model. Detailed definitions of variables can be found in reference 1 . 

In the framework of the fluid network approach, the entire fluid system is considered as consisting of 
a number of lumps connected by paths. Each lump has its own finite volume and is joined by path(s) to 
other lump(s), thus a network is formed. The fluid properties, such as density, temperature, pressure, 
viscosity, and mass fractions are computed and stored at the centers of lumps, while a unidirectional 
velocity is computed and stored in each path; its value can be less than zero, equal to zero or larger than 
zero. In general, the number of paths connecting to a lump can vary from lump to lump. In this work, the 
paths enter a lump are called inward paths, while the paths leave a lump are called outward paths. This 
definition facilitates the implementation of upwind differencing scheme and implicit time integration. 

In the following, the subscript for a path variable will be either j or i, and the subscript for a lump 

property will be i. This is the staggered-grid arrangement commonly used in many CFD schemes. In 
addition, a predictor-corrector algorithm is adopted to advance the temporal evolution of the variables. 

In the predictor step, a previously computed pressure is employed to solve the momentum equation 
for the path velocity, 

m it* - m jii'j . 

At + max ( w j > °X M j - U juw ) + 

max (- m j , 0 )(u* - u* dw ) = (p uwj - p dwj )Aj -K j] m j \ m j \ Aj + S uj 

( 6 ) 

where the subscript juw points to the upwind path of path j, while the subscript jdw points to the 
downwind path of path j. Conversely the subscript uwj points to the upwind lump of path j, while the 

subscript dwj points to the downwind lump of path j. Here, u, is the velocity of path j, mj is the mass 

flow rate of path j. Kg accounts for the laminar and turbulent viscous effects associated with the path j, it 
is a function of the geometry entities and the characteristics of the paths. For typical pipes, it is usually a 
function of Darcy friction factor. A comprehensive definition of Kjj is given in reference 6. Aj is the cross 
sectional area of path j. S UJ represents a generic source term such as a pump or additional forces exerting 
on a path. 
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Equation (6) can be recast into 


m i 


At 


■ K jj pj | m j | Aj + max ( m j , 0) + max(- m j , 0) 


n 

m i ll i 


■ max(mj,0)Uj UW + max(-mj,0)Ujd w + (p UW j-Pd W j)Aj + —^~ + S uj 


(V) 


In the subsequent discussion, we will use the variable Mg which is defined as 


• fYl j • • • 

Mtj - — - + KgPj | rtij | Aj + max (m y , 0) + max (- m y , 0) 
At JJ J J 

The energy equation is formulated in terms of the conservation of the enthalpy. 


m lv -m-h- 1 outward . inward 

11 1 1 ■+ J max (- mi , 0)(h* - h*j w £ ) - J max(m j ,0)(h* - h* wj ) 


At 


1 = 1 


7=1 


outward 

- I 

7=1 

inward 

* z 

7=1 


( ~Pdwi') + K fi m ^ + K H~r^'dwi ~ T i) + H f A wA T wall ~ T aw)i 

mi ^ 

^ max(mj , 0 ) Uj Aj )[(p dwj - Puwj ) + Kjj mj )] + Kj -j-(T uwj - 7] ) + HjA wj (T wall -T aw )j 

mj 


+S hi 


( 8 ) 


/* —2 / 1 

where the convective heat transfer coefficient is H = pCp u St , and 57 = — Pr is the Stanton 


number given by modified Reynolds analogy for turbulent flow in pipes, where /is the Fanning friction 
factor of the path, and Pr is the Prandtl number. T wa u is the specified wall temperature and T aw is the 
adiabatic wall temperature. L is the normal distance of a path. A H;7 is the wetted surface area of path j. S/j 
represents a generic source term for the energy equation such as the heat transfer from the surrounding 
solid. 

In the second, i.e., the corrector step, the continuity equation is formulated in terms of a pressure 
correction. It is assumed that the velocity is a function of the pressure gradient which, in turn, is a 
function of the two pressures on each side of a velocity path. Now consider a path j, which is connected to 
a upwind lump denoted by upj and a downwind lump denoted by dwj, applying a Taylor series expansion 
along the pressure coordinates, one has 


* dw ■ 
u i =u j + - 

"" uwj 


dp 


C Puwj-P i 


UWJ 


dll : 

) + ■ J 


d Pd 


C Pdwj Pdwj ) 


wj 


Using the definition of Mtj and equation (7), one finds 
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diij Aj diij Aj 

d Puwj - M t ) ’ M,J 

thus, the new velocity can be expressed as a function of the previously predicted velocity and the pressure 
change via 


u j= u j + 


A j 


( A Puwj - A Pdwi) 


WJ' 


M t j 


(9) 


Now applying the continuity equation to a lump yields 

outward inward 

\ 1 + Z Pe u t A e- Z p i ll i A i = 0 ( 10 ) 

£=1 ./=i 

Like other pressure-based methodology, the pressure change and the density change are related through a 
derivative evaluated at the condition of constant temperature: 


p-p + Ap = p + A p 


dp 

dp 


(ID 


where, for the perfect gas, 


dp 

dp 


J T 


- Py , and, for the high pressure fluid model, 


dp 

dp 


= 7_ 
J T a 2 


Substituting equations (9) and (11) into equation (10), and then neglecting the second order changes 
while applying the first- order-upwind differencing to the mass flux in the paths, one obtains the final form 
of the pressure correction equation 


V i Ap i 


outward 


R i T i At 7 =1 


m ax(-M, 0 ) + M 


( A Pi~ A Pdwl ) 


A * * a2 

Am j PjAi 

max(-L=L,0) + -L-L 

R j T j M t j 


M t l 
( A Pi~ A Puwj ) 


I 

1= 

inward 

+ 1 

j = 1 

outward r 

= Z max(-^M*,0)p^-max(^M*,0)p- 

e = i 

inward r 

+ Yj [ max ( A jUp°)Puwj ~ max(-Ajii*,Q)p* 
7=1 

m* - m- 
At 


(12) 


where V, is the volume of lump i. Once the pressure is obtained, the field values of the path velocity are 
corrected according to 
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( 13 ) 


n+ \ _ * : A (A Puwj A Pciwj ) 
ll j = Uj + Aj - 


and the enthalpy values are corrected according to 


h> r+l = /,*+— (XX 


outward 


IVVLl f LI - 

X 0 . 5 ( Api + Apd w ( ) A i ii] 


1 = 1 

inward 


vvuru i- 

X + A/7 mi , / )d / My +l ) 

j = 1 


where is the average of surrounding ^ at lump i. 

For each species s, there is one transport equation to be solved respectively 


yv) . — ffl-v ■ ou twa rd 

m iys,i m iy s ,i . V 1 / n+1 w+1 ^ 

7 -+ Zj m a x(-ms,0)(y ; -y s j wi ) 

A? U 

inward 

- X max(my,0)(y"f 
7=1 

X o r ( ys,dwi y s ,i ) 

^=i L ^ 

(ju + /j t ) A j 

X ^ (ys,uwj — y s ,i ) 

i=i L j 

+ ' F + A F 
1 &s,i 1 
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The fluid network based GDE is written as 


V i N kt~ V i N ki ou ^ rd ,+ 1 „+i 

Ar 1 = 1 


inward 


£ max( , O)(jV0 1 - jV^- ) 
7=1 


outward 

I 

7=1 

inward 

■1- I 


41. 


d,£ 


D, 


N 


M 


k,dwt 


4K 


d,j 


7=i L D> y 


N, 


k,uwj 


+v ; 


1 At 

i micro 

k = M b 


(16) 


where A* is the number density of particles in size bin k (the unit is #/m 3 ), k loops from 1 to the maximum 
number ( M B ) of sections for each particle size distribution, V d is the deposition velocity in a path, and D h 
the hydraulic diameter of a path. 


AN k 


At 


\micro 


is the summation of microphysical processes which contribute to the production and loss of the 


particles, 


AA h 

At 


micro 


AN, 


k,i 


At 


+ 


thermolphoresis 


AA hi 

At 


\coagidation 


+ 


AA ^ 
At 


+ 

nucleation 


A^V k ,i 


At 


condJevap . 


(17) 


Finally, the solution of all the transport equations leads to 

[Am = [B] (18) 

where [A] is a diagonal-dominated matrix, [. B ] is the source vector, and ()) is the solution vector. 
Equation (18) is solved here iteratively using a point-implicit relaxation method due to the unstructured 
nature of the fluid network. 


Experiment 

Recently, laboratory tests have been conducted to better understand the sources of the significant 
variations in the measured size and concentration of particles in the jet engine exhaust among several 
sampling and measurement approaches. The APEX-3 sampling line efficiency results are presented in 
reference 7. Line loss measurements for different tubing sizes and sample flow rates are reported in 
reference 8. 

In reference 7, the particle loss characteristics of the sampling lines used in the APEX-3 field 
measurements were examined in a laboratory environment. The goal was to quantify the particle 
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penetration efficiency as a function of the particle diameter. In these tests, monodisperse particles were 
generated by a tube furnace, and two well characterized particle instruments were used, one was 
connected to the entrance of the tested sampling line, the other to the exit of the tested line. The particle 
penetration efficiency of the sampling line is defined as the ratio of the measured downstream particle 
number density to the measured upstream particle number density. The layout of the test sampling lines 
has the same geometries as used in the field measurements. Two sets of data were obtained in the 
laboratory tests, one is for the case of the so called lm-probe line, i.e., the aerosol probe is located at 
1 meter downstream from the exit plane of the aircraft engine, the other one is for the case of the 30-m 
probe line. Schematics of the test setup can be found in the Appendix. 

In reference 8, the monodisperse particles were generated from a premixed flat flame burner. The 
particle penetration efficiency as a function of the particle diameter has been investigated for different 
material, length, diameter, and flow rate of the sampling line. A comprehensive line loss data sets are 
being compiled, and we have selected the test number 26 and test number 34 for modeling and 
simulation. The sampling lines are made of stainless steel, but they have different sizes, lengths, and flow 
rates. 


Assessment Results 

Four line loss experimental data sets are used for assessing the microphysical models currently in the 
NCC. These data sets are denoted as the APEX3 1-m probe sampling line (ref. 7), the APEX3 30-m probe 
sampling line (ref. 7), the UTRC test 26 sampling line (ref. 8), and the UTRC test 34 sampling line 
(ref. 8). Monodisperse particles have been used in all of these tests. And, for the flow conditions and the 
particle size range tested, the experimental data strongly suggest that diffusion of particles to the walls is 
the primary particle loss mechanism, other effects, such as the bends in sampling line, are negligible. 

Three different levels of simulation have been carried out: CFD-based simulation, network-based 
simulation, and lower-order temporal simulation. It should be noted here that all of them use the same 
chemical and microphysical models. Each calculation requires the specification of particle size 
distribution and concentration at the entrance of the sampling line, and, instead of monodisperse particles, 
we have used polydisperse particles for all of the numerical simulations. The 12-bin discrete distribution 
of these polydisperse particles is given in table 1 . Calculations in which all of the chemistry and 
microphysical models for the particles were “turned on” have been performed, and the results are very 
much the same as those from the calculations which include only the effects of particle diffusion. This is 
consistent with the experimental suggestion that diffusion of particles to the walls is the primary loss 
mechanism in the cases considered, and other effects, such as coagulation among particles, are 
insignificant. It is mentioned here that the results presented in the following are obtained by using 
polydisperse particles with the option of coagulation being “turned off’ in the calculations. 


TABLE 1.— DISTRIBUTION OF PARTICLE SIZE AND NUMBER DENSITY 


Bin number 

Diameter, 

Number density. 


(nm) 

(#/m 3 ) 

1 

3 

8585090 

2 

4.48 

453587191 

3 

6.7 

10742262056 

4 

10 

95160588390 

5 

15 

315315201902 

6 

22.4 

390804103104 

7 

33.5 

181175611146 

8 

50 

31417148866 

9 

74.8 

2037793720 

10 

112 

49440184 

11 

167 

448669 

12 

250 

1220 
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1-m particle probe sampling line 

Since it is known that the effect of bends in the line on the particle loss is insignificant in this case, we 
consider a straight pipe consisting of two segments with different diameters. The total length of the 
sampling line is 23.5 m. The internal diameter of the first 9.5 m segment is 0.009525 m (3/8 in.). The 
internal diameter of the second 14 m segment is 0.019050m (3/4 in.). The CFD-based simulation assumes 
the flow is axisymmetric, and a total of 57,680 quadrilaterals are used to discretize the computational 
domain. Figure 1 shows a part of the mesh. The grid distribution used to conduct the lower-order 
temporal simulation is shown in figure 2, the number of grid points is 40. Similar arrangement is used in 
the fluid network simulation, i.e., there are 40 lumps and 39 paths. 

The incoming flow rate is 50 liter per minute; hence, the equivalent inlet velocity is 1 1.7 m/s. The 
temperature of the carrier gas is 295 K. The upstream pressure is 740 Torr (98627.2 Pa). The downstream 
pressure is 690 Torr (91963.2 Pa). Computed results and the experimental data in terms of the particle 
penetration efficiency are presented in figure 3. The parameter “a” used in the CFD-based simulation 
accounts for the sticking effect of the wall. The value of a = 0.5 means that the number density of the 
particles at the wall is half of the value at the nearest cell. Figure 3 indicates that, for smaller particles 
(say, d< 10-20 nm), temporal and fluid network methods over-predict the penetration efficiency, while 
CFD-based method under-predicts the penetration efficiency. For larger particles, all three methods 
provide reasonably good results. 



Figure 1 . — Schematic of the mesh used in the CFD-based simulation. 


Q Q Q Q Q O □ Q G G G G G G G G G G G G G G G G G G G G O 


Figure 2. — Schematic of the grid distribution used in the lower-order temporal simulation 


APEX3 1 m-probe line loss (Monodisperse aerosol source) 

T=295K. Pupstream=740 Torr, Pdownstream=690 Torr. Flow=50 lpm( 1 1 .7 m/s) 



Figure 3. — Line loss comparison: flow rate = 50 LPM, inside diameter = 0.009525 m for 
the first 9.5 m segment, inside diameter = 0.01905 m for the second 14 m segment, 
T g as=295K, Pout = 91963 Pa. 
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30-m particle probe sampling line 

The simulated sampling line is a straight pipe with an internal diameter of 0.01905 m (3/4 in.). The 
total length of the sampling line is 35 m. The incoming flow rate is 50 liter per minute. Thus, the 
equivalent inlet velocity is 2.9 m/s. The temperature of the carrier gas is 295 K. The upstream pressure is 
740 Torr (98627.2 Pa), and the downstream pressure is 720 Torr (95961.6 Pa). Computed penetration 
efficiency of the particles and its experimental counterpart are presented in figure 4. For smaller particles 
(say, r/<30-40 nm), all methods, except for the fluid network method with modified slip factor, over- 
predict the penetration efficiency. For larger particles, all methods yield reasonably good results. 

The number of grid points used in the lower-order temporal simulation is 50. The fluid network 
simulation uses 50 lumps and 49 paths. Two different meshes are used for the CFD-based simulation: 
75,040 quadrilaterals used in an axisymmetric simulation and 380,160 hexahedrons used in a three- 
dimensional simulation. 

In the CFD-based simulations, in stead of using the gradient of particle number density to compute 
the molecular diffusion near the wall, we now take the approach of using the deposition velocity to 
compute the diffusion flux near the wall. Consequently, the value of the parameter “a” is now set to zero. 

Figure 4 also demonstrates the effect of the so called slip factor on the results obtained from the fluid 
network simulation. Typically, the Cunningham-slip factor is defined as 

Ccj =1 + &z ; (1.257 + 0.4e~ 1,VKni ) (19) 

where Kn , is the Knudsen number defined as the ratio of the mean free path of the carrier gas, /.„, to the 
diameter of the size-/ particle 


APEX3 30m-probe line loss (Monodisperse aerosol source) 

T=295K, Pupstream=74Q Torr, Pdownstream=720 Torr, Flow-50 lpm(2.9 m/s) 



Figure 4. — Line loss comparison: flow rate = 50 LPM, inside diameter = 0.019050m, 
Length = 35 m, T gas = 299K, P ou t = 95961 Pa. 
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( 20 ) 


Krii = 



The mean free path of the carrier gas depends upon the pressure, temperature, molecular weight of the 
carrier gas, and the universal gas constant. 

It is observed that, for smaller particles, equation (19) leads to values of Cc, which are much larger 
than one, while for larger particles, the values are close to one. To allow more control of the values of this 
empirical parameter, equation (19) is modified as 

Cc[ = 1 + Ktij (1 .257 + 0.4e -1 • : {/Kn > ) x b f (21) 

where b, is an particle-size dependent parameter. Experiences show that, generally speaking, when b, is 
larger than one, it leads to a reduction in the computed penetration efficiency. In the present simulation, 
values between 1.0 and 1.2 are specified for b h depending on the particle size range. 

UTRC Test 26 sampling line 

The tested sampling line is a straight pipe having an internal diameter of 0.008 1 026 m (0.3 1 9 in.). 

The total length of the sampling line is 15.24 m. The incoming flow rate is 50 liter per minute. Thus, the 
equivalent inlet velocity is 16.2 m/s. The temperature of the carrier gas is 298.6 K. The upstream pressure 
is 648.9 Torr (86512 Pa), and the downstream pressure is 540.2 Torr (72020 Pa). 

The CFD-based simulation assumes axisymmetric flow, and 34,720 quadrilaterals are used for the 
mesh. The number of grid points used in the lower-order temporal simulation is 3 1 , the fluid network 
simulation uses 31 lumps and 30 paths. Computed as well as measured penetration efficiency are shown 
in figure 5. For smaller particles (say, d<10 nm), all numerical methods predict similar results. For larger 
particles, both CFD-based simulation and lower-order temporal simulation under-predict the penetration 
efficiency with respect to the experimental data. The results from the fluid network simulations compare 
well with the measured data. 


UTRC Test 34 sampling line 

The tested sampling line is a straight pipe having an internal diameter of 0.01 1684 m (0.460 in.). The 
total length of the sampling line is 30.48 m. The incoming flow rate is 72.5 liter per minute; hence, the 
equivalent velocity is 1 1.3 m/s. The temperature of the carrier gas is 299.9 K. The upstream pressure is 
632.5 Torr (84324 Pa), and the downstream pressure is 551.6 Torr (73539 Pa). 

The CFD-based simulation is assumed to be axisymmetric, and the mesh has 75,040 quadrilaterals. 
The lower-order temporal simulation uses 46 grid points. The fluid network simulation uses 46 lumps and 
45 paths. Computed and measured penetration efficiencies are presented in figure 6. The predictions from 
the fluid network simulation with modified slip factor and the lower-order temporal simulation agree with 
the experimental data very well. The result from the fluid network simulation without modified slip factor 
over-predicts the penetration efficiency. 
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UTRC line loss (Monodisperse aerosol source ) 

T=298.6K. [-’upstream =64 8.9 Torr. Pdownstream=540.2 Torr, Flow=50 lpm( 16.2 m/s) 



Figure 5. — Line loss comparison: flow rate = 50 LPM, inside diameter = 0.00810265 m, Length 
1 5.24 m, T gas = 298.6 K, P 0 ut = 72042 Pa. 


UTRC line loss (Monodisperse aerosol source) 

T=299.9K. Pu pstream=6.32.5 Torr. P<Jownstream=551.6 Torr, Flow=72.5 lpm( 1 1.3 m/s) 



Figure 6. — Line loss comparison: flow rate = 72.5 LPM, inside diameter = 0.01 1684 m, Length = 
30.48 m, Tgas = 299.9 K, P ou t = 73539 Pa. 
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Concluding Remarks 

A fluid network simulation capability is now available in the NCC. Together with the previously 
available CFD-based simulation capability and the lower-order temporal simulation capability, we now 
have a suite of computational tools which are versatile and cost effective for the modeling and simulation 
of transport of aerosols in the various flow paths, such as the engine internal passages, the sampling 
probes, the sampling lines, and the plumes. The focus of the present effort is the assessment of these three 
levels of simulation when applied to the prediction of particle loss in the sampling lines used in the 
measurement of particulate matter size and concentration in high-temperature, high-velocity jet engine 
exhaust. The predicted particle penetration efficiencies have been compared with the corresponding 
experimental data from four different laboratory tests. Generally speaking, the predicted results agree well 
with the measured data, when the particle diameters are larger than, say, 30 nm. Based on the results 
obtained so far, among the three simulation approaches, the fluid network approach with modified slip 
factor seems to be the optimum choice for the modeling and simulation of the sampling line system. In 
addition to help establish standard methodologies for measuring the size and concentration of particles in 
the jet engine exhaust, the present work also represents the first step of a long term effort of validating 
physics-based tools for the prediction of aircraft particulate emissions. Obviously, this requires the 
acquirement of particulate emissions data specifically for the purpose of validating the microphysical 
models for aerosols. 
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Appendix — Laboratory Test of APEX-3 Sampling Lines 

To determine the particle loss characteristics of the sampling lines used in the APEX-3 field 
measurements, laboratory tests were conducted to quantify the particle penetration efficiency as a 
function of the particle diameter. Monodisperse particles were generated by a tube furnace. Two well 
characterized particle instruments are used, one is connected to the entrance of the tested sampling line, 
the other is connected to the exit of the tested line. The particle penetration efficiency of the sampling line 
is the ratio of the measured downstream particle number density to the measured upstream particle 
number density. 

The layout of the test sampling lines has the same geometries as used in the field measurements. Two 
sets of data were obtained in the laboratory tests, one is for the case of the lm-probe line, i.e., the aerosol 
probe is located at 1 m downstream from the exit plane of the aircraft engine, the other one is for the case 
of the 30-m probe line. 

The following schematics are taken from reference 7 to illustrate the laboratory line test setup, the 
layout and geometry of the lm-probe line, and the layout and geometry of the 30-m probe line: 

Laboratory test 


Line Loss Test Setup 


System was carefully 
leak-checked to 
minimize stray particle 
counts to allow going 
to high sample 
dilutions/flow rates 


Monodisperse Aerosol 
Source (NaCI from Tube 


Furnace) 


ej] — Dilution Gas 


Sampling line 
to be tested 



Particle instruments 
were tested side-by- 
side prior to 
experiment to 
determine any 
differences in size 
sensitivity or 
response. 


System flow was 
adjusted to simulate 
actual sampling 
conditions and to 
evaluate losses as a 
function of flow rate. 


Distribution Manifold / 

rA i 


Jn MFM 
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APEX-3 1-m probe line 


APEX-3, Aerosol Sampling Line 



APEX-3 30-m probe line 
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